Scale Invariant Fractal and Slow Dynamics in Nucleation and Growth Processes 
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We propose a stochastic counterpart of the classical Kolmogorov-Johnson-Mehl-Avrami (KJMA) 
model to describe the nucleation-and-growth phenomena of a stable phase (S-phase). We report 
that for growth velocity of S-phase v = s(t)/t where s(t) is the mean value of the interval size x 
of metastable phase (M-phase) and for v — x/t(x) where t(x) is the mean nucleation time, the 
system exhibits a power law decay of M-phase. We also find that the resulting structure exhibits 
self-similarity and can be best described as a fractal. Interestingly, the fractal dimension df helps 
generalising the exponent (1 + df) of the power-law decay. However, when either v — Vo (constant) 
or v = o~/t (a is a constant) the decay is exponential and it is accompanied by the violation of 
scaling. 
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PACS numbers: 68.55.Ac, 64.70.Kb, 61.50.Ks 



The kinetics of phase transformation via nucleation 
and growth of a stable phase (S-phase) represents one 
of the most fundamental topic of interest in both sci- 
ence and technology pj. It plays a key role in metal- 
lurgical applications as well as in many seemingly unre- 
lated fields of research. The kinetics of crystal growth 
, the domain switching phenomena in ferroelectrics Q 
and the DNA replication in organisms are just a few 
examples to name. Much of our theoretical understand- 
ing of such phenomenon is provided by the Kolmogorov- 
Johnson-Mehl-Avrami (KJMA) model. It has been for- 
mulated independently by Kolmogorov, by Johnson and 
Mehl and by Avrami in and around the 1940s 0- It still 
remains one of the most studied theory and is the only 
means of interpreting the experimental data to gain in- 
sight of the process. The basic algorithm of the KJMA 
theory is trivially simple. One time unit of the process 
in one dimension can be defined as follows. 

i) Randomly a position is selected for nucleation of 
point-like seed of S-phase on an interval of M-phase. 

ii) Upon nucleation the seed starts growing isotropi- 
cally on either side with a constant velocity vq at 
the expense of the decay of M-phase. 

Hi) Whenever two such growing phases from opposite 
sides meet, growth ceases immediately at the point 
of contact, while continuing elesewhere. 

iv) The process is repeated ad infinitum. 

In the context of nucleation and growth phenomena, 
the question we usually ask is: What is the fraction of 
M-phase that still survives at time tl According to the 
KJMA theory, such quantity in d dimensions obey the 
exponential decay known as the Kolmogorov- Avrami law 



$(<) = cxp 
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where, V decribes the constant nucleation rate per unit 
volume and f2<j is the constant volume factor of the d 



dimensional hypersphere (f2<j = l,7r,47r/3 for d — 1,2,3 
respectively and d+ 1 is often called as the Avrami expo- 
nent) Q. Recently, the derivation of correlation function 
and its connection to the scattering cross-section p , the- 
ory of grain-size populations Q , and the generalization of 
the KJMA theory to multiple stable phase has made a 
significant contribution as it now provide a better means 
of interpreting the experimental data. Alongside, there 
have been reports too that the experimental data in some 
cases do not fit a straight line in the plots of loglog[<j>(f)] 
against log[£] [9j,ll0(. It clearly raises some concerns and 
hence it requires further theoretical interest in order to 
find varients of the model which would be suitable under 
different physical situations. 

In this paper, we lift the restriction imposed by the rule 
(ii), the conventional constant growth velocity, by choos- 
ing various form of velocities. Precisely, we choose veloci- 
ties that decreases as the system proceeds and check how 
such decreasing grwoth velocities of different strength in- 
fluence the resulting dynamics. In fact, a more realistic 
theory does require the inclusion of a growth velocity 
that may change in the course of time. The idea to move 
away from the constant growth velocity is born by the 
following observation. The problem in question is intrin- 
sically a kinetic one and hence, as the time proceeds, the 
size and shape of both the phases change continuously 
altering the strength of their competition between the 
surface tension and the free energy density difference be- 
tween the two phases. The strength of this competition 
is the one that ultimately determines the growth velocity. 
This implies that the growth velocity of S-phase should 
therefore depend either on time or on the size and shape 
of one, or both, of the phases or on proper combinations 
of all these parameters in some sense. The inlcusion of 
velocity, that is not constant, appears to reveal a num- 
ber of new interesting results. The most significant of all 
is the emergence of the scale invariant fractal geometry 
accompanied by a power-law decay of M-phase instead 
of the traditional exponential decay. We show that the 
system exhibits such fascinating behaviour only if the 
growth velocity decreases as the system eveolves. To be 
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more precise, the systems exhibit scale-free fractal if the 
growth velocity at any given time t depends on the mean 
size s(t) of M-phase i.e., v = s(t)/t or if it depends on 
the mean nucleation time t(x) i.e., v = x/t{x). 

In this work, we mainly focus on extracting the key 
features of the phenomena (at least thier qualitative be- 
haviour) with the least possible mathematical complex- 
ity. We therefore choose the one-dimensional (Id) system 
as we can solve it for some special cases analytically and 
exactly. The Id model can in fact fully capture the qual- 
itative bahaviour of the key quantities of interest namely 
the decay of M-phase. The Kolmogorov-Avrami formula 
itself is a testament to its justification which clearly shows 
that the exponential decay of M-phase is common to all 
dimensions. The distribution function C(x, t) of M-phase 
of size x at time t for random nucleation and growth 
mechanism evolves according to the following integro- 
differential equation 

dC(x,t) 
dt 



= -xC(x,t)+2 / C{y,t)dy (2) 



d_ 

dx 



'v(x,t)C(x,t)]. 



The first two terms on the right hand side represent the 
decrease and increase of interval size x of M-phase due 
to the nucleation of S-phase on x and y > x respectively. 
Every nucleation event creates two new intervals. The 
factor '2' in the integral term implies that either of the 
two new intervals can be of size x. Finally, the third term 
on the right hand side is derived by taking into account 
all the possible ways the distribution function C(x, t) still 
remains in the size range [x, x + dx] despite loss and gain 
of x due to continuous growth in a span of infinitesimal 
time dt. For instance, if we assume that the process is 
homogeneous in time then all the possibilities are 

COM + ^)l gr owth = [l-v(x,t)dt/dx]C(x,t)+ (3) 

[v(x + dx, t)dt/dx]C(x + dx, t), 

which is in fact the third term on the right of Eq. J2J). 
Here, v{x,t)dt/ dx (v(x + dx,t)dt/ dx) is the fraction of 
concentration C(x, t) that is lost (gained) in time dt due 
to the growth of S-phase. 

We now attempt to solve Eq. J5J) subject to the initial 
conditions 



C(x,0) = -5(x 

Li 



lim 



C(x 1 0)dx = 0. (4) 



This ensures that there is no seed of S-phase at t = 0. 
Once we know C(x, t), we can immediately find the frac- 
tion of M-phase that remained untransformed <£>(£) = 
Jp 00 xC(x, t)dx and the number density of the intervals of 
M-phase N(t) — C(x, t)dx. The fraction of M-phase 
covered by S-phase is related to $(t) via 6{t) = 1 — $(t) 
that evolves as 



dO(t) 
dt 



v(x, t)C(x, t)dx. 



(5) 



The quantities $(£) and N(t) can also be used to know 
how the average interval size s(t) = &(t)/N(t) varies in 
time. 

The two parameters involved in the definition of ve- 
locity are length and time. The growth velocities that 
we can choose are as follows: (i) v(x,t) — vq, the tradi- 
tional constant growth velocity (classical KJMA model), 
(ii) v(x,t) = a/t where a is a constant and bears the 
dimension of length, (iii) v(x,t) = s(t)/t where s(t) is 
the mean interval size of M-phase at time t and finally 
(iv) v(x, i) — mx/r(x) where t(x) is the mean nucleation 
time and to is a dimensionless positive constant 11]. 

For growth velocities (a) - (c), which are independent 
of x, the solution of Eq. J5J) can be obtained by substi- 
tuting the ansatz 



C(x,t) = B(t)exp[-xt]. 



(6) 



Here the time dependent prefactor B(t) obeys the ordi- 
nary differential equation 



d\nB{t) 
dt 



2/t - v(t)t. 



(7) 



We have to solve Eq. JJJ) subject to the initial condition 
B(0) = 0, followed by Eq. (0J. To proceed further and 
for clarity, we will treat each case independently. 

(a) Solving Eq. J7J for v = vq we obtain B(t) — 
t 2 exp[— vot 2 /2] and therefore 



C(x, t) = t 2 exp[-xt - v t 2 /2]. 



(8) 



Using this in the definition of <&(£), we immediately find 
$(t) = exp[— vot 2 /2] which is the celebrated Kolmogorov- 
Avrami formula in one dimension, except the factor T. 
Note that Eq. J2J describes the random sequential nu- 
cleation of one seed at each time step and hence T = 1. 
It clearly demonstrates not only how well the theory 
works but also its simplicity. Furthermore, we find that 
the number density of intervals of M-phase changes as 
N(t) - te~ v ° t2/2 . It reveals that N(t) rises linearly at 
the early stage due to the nucleation of S-phase. At the 
late stage on the other hand, N(t) decreases exponen- 
tially due to the fast coalescence of the neighbouring sta- 
ble phases 01 ■ The mean interval size s(t) on the other 
hand decays as s(t) ~ t~ x . 

(b) We now solve Eq. (J7J) for the velocity v(t) — a/t. 
to give 



C(x, t) = t 2 exp[-(z + a)i\. 



0) 



In this case too, we find that the fraction of M-phase de- 
cays exponentailly $(i) = exp[— at] but slower than that 
for constant velocity. Here the Avrami exponent a = 1 
which corresponds to the one typically known for the het- 
erogeneous nucleation and growth processes despite 
the fact in the present case it strictly describes the homo- 
geneous nucleation. It proves that the Avrami exponent 
not only depends on the nature of the nucleation process 
but it also depends on the nature of the growth velocity. 
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We also find that the number density of M-phase varies 
as N(t) ~ te~ at . This reveals that at the early stage, 
it rises linearly like the model (a). However, at the late 
stage, the coalescence takes place less frequently than in 
the model (a). 

(c) In order to be able to solve Eq. for v(x,t) = 
s(t)/t, we first need to know explicitly how the mean in- 
terval size s(t) decays in time. We have already solved 
Eq. (J2J for two different growth velocities and in both 
cases we found that the mean interval s(t) decreases with 
time as s(t) ~ t~ l . Interestingly, it is also the case with- 
out the growth term which is known as random scission 
model [lj|. All this suggests that s(i) = kt" 1 is the most 
generic choice for the mean interval size, where k is a 
dimensionless parameter. Incorporating it into the defin- 
tion (c) gives v(x, t) — k/t 2 . Solving Eq. J7J for this 
growth velocity yields 



C(>,*) = t 



2-k 



exp[ 



rt}. 



(10) 



Unlike the previous two cases, it is of particular interest 
for the following reasons. We find that all the moments, 
M n (t) of C(x,t) where M n (t) = f™ x n C(x,t)dx, exhibit 
a power-law behaviour M n (t) ~ f _ ( n_ ( 1_fc )) and there- 
fore, 



$(i) = t- k . 



(11) 



The emergence of such a slow decay makes it a good 
candidate to look for a scale invariance in the system. 
The k value here plays the crucial role in determining 
the dynamics of the system and is bound by the lower 
and the upper limits and they are fixed by the following 
physical constraints. The lower bound is fixed by the 
behaviour of $(t) that should be increasing function of 
time and hence it demands k > 0. The upper bound on 
the other hand is fixed by the constraint that the number 
density N ~ t x ~ k should be an increasing function of 
time, at least, during the early and intermediate stage. 
This immediately provides the upper bound k < 1. The 
only non-trivial and physically interesting k value are the 
ones that stay within the interval [0, 1]. In passing, we 
note that in this model too the average interval size s(t) 
decays as s(t) = t _1 . 

(d) Finally, we consider the case where v(x,t) = 
mx/r^x), however, according to Eq.Q, the typical time 
r between two nucleation events on an interval of size x 
is t = x^ 1 and therefore v(x,t) = mx 2 revealing that 
the growth velocity decreases increasingly fast as the size 
of M-phase decreases since x itself is a decreasing quan- 
tity with time. For simplicity reason, we only consider 
the case m = 1. Invoking a simple dimensional analysis 
allows us the following scaling ansatz 



C(x,t) ~t w (f>(xt), 



(12) 



such that t w bears the dimension of C(x,t) and 0(0 is 
the scaling function. Upon substituting the ansatz into 
Eq. an d seeking the self-similar solution immediately 



reduces the original partial integro-differential equation 
into an ordinary hypergeometric differential equation for 



«i-0^ 



{(l + «0-3£}^p + 0(0=0. (13) 



The most nontrivial feature of 0(0 is its universality, 
in the sense that it is independent of the initial con- 
dition. The only physically acceptable solution of the 
above equation is 

0(0 = 2 Fi(l + >/2,l- V2; l + «;;0, (14) 

and asymptotically it behaves as 0(0 ~ exp[— (V2 — 1)£] 
|15|. The yet undetermined exponent w of the scaling 
ansatz is fixed by the conservation principle. However, as 
the intervals of M-phase are lost to S-phase, the system 
violates the simple mass conservation. Nevertheless, it 
still obeys a conservation law although this time it is a 
non-trivial quantity as we shall see below. To do so, we 
incorporate the definition of M n (t) into Eq. J5J and then, 
after some simple algebraic manipulations, we are able to 
write the following rate equation for M n (t) 



dM n (t) {n-a)(n + a + 2) 
Jt "~ (n+ 1) 



M n+1 (t). 



(15) 



Thereafter, demanding the steady state condition 
dM n (t)/dt — immediately reveals that it is the value 
a = (V2 — 1) that fixes the exponent w = \/2. The 
asymptotic solution for C(x,t) is therefore, 



COM) 



t^expf 



-(V2-l)xt\. 



(16) 



We thus find that <&(i) decays algebraically with a non- 
trivial exponent 0.5857 whereas N(t) grows with an ex- 
ponent 0.4142. However, the average interval of M- 
phase s(t) decays exactly like all previous three cases i.e., 
s(t) ~ t~ x . For model (c) too, we could substitute the 
ansatz Eq. 112|l into Eq. J5J) and find that the scaling 
function 0(0 obeys an ordinary differential equation. Its 
solution would have been the same as we obtained by us- 
ing the direct method, Eq. (fTU|) . This further supports 
the fact that the model (c) admits scaling. However, 
applying the same treatment to the model (a) and (b) 
would immediately reveal that neither of these two cases 
yield a reduction to an ordinary differential equation of 
one variable. This suggests that the model (a) and (b) 
do not admit a scaling solution. In fact, the exponential 
decay of M-phase, is a clear signature of the violation 
of scaling. It is note worthy to mention that the number 
density in (c) and (d) increases for all time - a sharp con- 
trast with the models (a) and (b) where N(i) increases 
only at the early stage. This is due to the fact that in 
models (c) and (d), the growth velocity decreases in time 
in such a way that two growing phases from opposite di- 
rection hardly coalesce. In fact, the growth of S-phase 
virtually stops prior to coalescence. 
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The existence of scaling and the non-trivial exponent 
of $(£) provides an extra motivation to go beyond the 
simple scaling description. This can be done by invoking 
the idea of fractal analysis, as it has been a very useful 
tool to obtain a global exponent called fractal dimension. 
To do so, we need a proper yardstick to measure the 
size of the set created in the long time limit. The most 
convenient one is the average interval size s(t) ~ t . 
Using this yardstick, we find that the number of s(t) we 
need to cover the system scales as 

lim N(s(t))~s(t)- d f. (17) 

s(t) >0 

The exponent df is known as the fractal dimension or 
the Hausdorff-Besicovitch dimension ^(| of the resulting 
set created by the nucleation and growth process. For 
v(x, t) = k/t 2 , we find df = (1 — k) and v (x, t) = x 2 yields 
df = 0.4142. For general m in (d), the fractal dimension 
is d f (m) = + 1/ni) + ± ^/(l + 1/m) 2 + 4/m (for 

further details see [13). In both cases ((c) and (d)), 
the Kolmogorov-Avrami formula can be replaced by the 
following general power-law decay of M-phase 

~ i-d-<*/). (18) 

It reveals a generalised exponent (1 — df) that can quan- 
tify the extent of decay of M-phase. The coverage 9(t) by 
the S-phase reaches its asymptotic value 6>(oo) = 1 again 
follow a power-law 

9(oo)-9(t) ~t-( 1 ~ d f\ (19) 

which is reminscent of the Feder's law in RSA [Tsj . 

To put it into a context, we have studied the random 
nucleation and growth processes of a stable phase for dif- 
ferent growth velocities. So far, the KJMA model has 
only been been studied with constant velocity. However, 



we conceive in systems where the dynamics is governed 
by the nucleation and growth processes, it is more logi- 
cal to consider that the growth velocity of S-phase may 
depend on the size of M-phase present at every instant of 
time. Invoking the idea of such velocites, we have found 
the following interesting and previously unknown results. 
The average interval size s(t) of M-phase shows a power- 
law decay s(t) ~ t^ 1 irrespective of the detailed choice 
of the velocity. The fraction of M-phase that still survive 
at any given time is very sensitive to the specific choice 
of the velocity. The resulting structure exhibits a scale- 
invariant fractal if only the growth velocity depends on 
either the mean size of the M-phase s(t) or on the mean 
nucleation time t(x) at every instant of time. Moreover, 
only in such cases the system show a power-law decay of 
the M-phase and we hve shown that exponent of the de- 
cay law can be generalised (1 — df) by incorporating the 
fractal dimension df. However, in the case where both 
or one of the parameter of the definition of velocity is as- 
sumed a constant value (e.g., v = Vq or v = constant/i) 
the decay is exponential in nature which is always ac- 
companied by the violation of scaling. 

We conclude with the following words. The power-law 
decay of M-phase can be assumed to be a generalized 
formula replacing the classical Kolmogorov-Avrami law, 
provided, the distribution of M-phase in the late stage 
describes a scale- free fractal. The fractal diemsnion df is 
the quantitative measure of the notion that the density of 
M-phase is less at larger length scale. We believe that the 
present work will have a significant impact in changing 
the way we intended to interpret the experimental data 
as we are now aware of the fascinating results due to the 
decelerating growth velocity. 
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